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O ; Abstract 

Forward-backward correlations are explored within the two-component clan model 
of multiparticle production. It is found that existing data are well described, and, in 
hh collisions, that clans must be allowed to leak particles from one hemisphere to 
the other. General formulae given for the symmetric case are then extended to the 
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asymmetric one, which is relevant for pA and AB collisions. 



The analysis of multiplicity correlations between two different regions of phase-space pro- 
vides important information on the dynamics of multiparticle production processes. 

Due to space limitations, the reader is referred for all details and references to papers [|l]j 
and®. 

I will at first consider symmetric reactions (like e + e~ annihilation, pp and pp collisions) 
with a symmetric choice of the forward (F) and backward (B) regions: in e + e~ annihilations 
the F regions is chosen randomly between those defined by a plane through the collision point 
and perpendicular to the thrust axis; in hh collisions the collision axis is chosen as reference. 
B is always the symmetric region. The forward-backward (FB) correlation strength, b FB , is 
defined as usual: 

((n F - n F )(n B - n B )) 



'FB 



(1) 



[((n F -fi F y)((n B -n B y)] 1/a ' 

where n F if the number of particles in F and n B the number of particles in B. A more precise 
study of the correlations investigates the dependence of the average number of particles in one 
hemisphere on the number of particles in the opposite one, say n F (n B ). 
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In hh collisions the correlation strength grows logarithmically with the cm. energy ^fs and 
is rather large (from 0.16 at 63 GeV to 0.7 at 1800 GeV); in e + e _ annihilation the correlation 
strength is rather small (0.1 at 91 GeV). Furthermore, b FB was measured in 2-jet and 3-jet 
events separately at LEP1 and was found to be compatible with zero in both cases. 

The mechanism of the weighted superposition of different classes of events, which was 
found to describe the shoulder structure of multiplicity distributions and the oscillations of their 
H q moments, can be tested on FB correlations. The classes of events are defined according 
to the reaction considered: in hh collisions we define soft events (events without mini-jets) as 
class 1, and semi-hard events (events with mini-jets) as class 2; in e + e~ annihilation, class 1 are 
two-jet events, class 2 are three-jet events. 

The starting point is the weighted superposition equation for the joint distribution of n F and 
n B particles: 

P(n F , n B ) = aPi(n F , n B ) + (1 - a)P 2 (n F , n B ), (2) 

where a is the fraction of events of class 1. Based on the above equation and symmetry proper- 
ties, it is possible to obtain explicitly a general formula for b FB , which depends of course on the 
first two moments of the general multiplicity distribution (MD) in each class i = 1,2 (namely 
the average multiplicity fii and the dispersion D n {) and on the correlation strength in each class, 
bi. 

ahDl^l + b 2 ) + (1 - a)b 2 D 2 n2 {l + bj) + \a(l - a){n 2 - ntfjl + + b 2 ) 
FB ~ aDl x (l + b 2 ) + (1 - a)D 2 n 2 (1 + b x ) + ±a(l - a)(n 2 - ^i) 2 (l + &i)(l + b 2 ) 

(3) 

We apply first our result to OPAL data, setting b\ = b 2 = as the data suggest. We take 
the values of the moments of the MD in each class from previous fits and obtain b FB = 0.101 
to be compared with the OPAL result 0.103 ± 0.007. Thus we conclude that the weighted 
superposition mechanism describes well FB correlations in e + e~ annihilation. 

We now turn to pp collisions at 546 GeV (UA5); the experimental result is 0.58 ± 0.01. 
Taking the shoulder description from a fit to the overall MD, one finds that Eq. (Q) with bi = 
b 2 = gives a too small result, indicating that FB correlations are needed in each substructure. 
On the other hand, if particles in each substructure were totally uncorrelated, one would get a 
result which is much larger than experimentally found. We conclude that particles are emitted 
with partial correlation, as already suggested by UA5. 

Since we have succesfully used a two-step production process to describe multiparticle dy- 
namics, leading to the negative binomial (NB) shape of the MD, and which naturally accom- 
modates partial correlation, we apply it to the study of FB correlations in a generalised way: 

- In the first step, N independent objects (generalised clans) are produced according to a 
given MD, V{N). 

- In the second step, particles are produced within each clan according to a MD Q(n c ). 

This two-step process is used to describe each component (with possibly different parameters 
in each step): from now on, all our formulae refer to each component separately, but we drop 
the component index in the formulae. We further assume that clans are produced independently 
in rapidity thus the distribution in the number of forward clans, N F , at fixed total number of 
clan, N, is binomial. The successful NB parametrisation of the data is obtained when V(N) is 
the Poisson distribution and Q(n c ) is the logarithmic distribution. 
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Consider now the situation in which particles produced by an F (B) clan all remain in the F 
(B) region: if clans follow a Poisson MD, as is implied by the very good fit to the data obtained 
with the NB distribution (NBD), one obtains b = 0, independently of the details of the MD 
inside clans. It is therefore necessary to allow clans to leak particles from one hemisphere to 
the other. In order to do this, a new parameter is introduced, p, which is defined as the average 
fraction of particles in a clan which remain in the same region and do not leak away to the 
opposite one. As clans are classified F or B according to where the majority of the particles fall, 
one has 0.5 <p<l. We also define q — 1 — p. 

The general result is now 

b = D 2 Jn - D 2 Jn c - (p - gfn c + 4 7 /n c 
Dl/n + Dl/n c + (p - q) 2 n c - A^/n c ' 

where 7 is the covariance between F and B multiplicities within a clan. The above formula 
simplifies considerably in the case of the NBD, and even further if we assume that particles 
within a clan fall in F or B independently, since then 7 = (D 2 — n c )pq: 

' Wn (5) 



1 - 2ppq 



where (3 is related to the average number of particles per clan: n c = (3 [((3 — 1) ln(l — /3)]~ . 

We can now use experimental data to extract the values of the leakage parameter p which 
we do not know, then we use the NBD structure to reproduce the almost linear relation between 
n F and n B . Since at 63 GeV the mini-jet sample can be considered negligible, we use only 
one component and obtain p soit = 0.78. Because the number of particles per clans in the soft 
component varies very little from 63 to 900 GeV, we assume now that also the leakage parameter 
is essentially constant in the GeV region. Using this value for p soft and the measured value of 
the correlation strength at 900 GeV, we can now calculate p S emi-hard = 0.77. The just mentioned 
values can then be used to predict the behaviour of the correlation strength with energy: e.g., it 
can be shown that it will not continue to grow linearly with In */s in the multi-TeV region if the 
leakage parameter J9 sem i-hard remains constant [[[]] . 

Within the same two-step process framework, we can extend our formulae to the case in 
which either the reaction is asymmetric (e.g. consider proton-nucleus or nucleus A-nucleus B 
with A 7^ B) or the F and B regions are asymmetrically defined, or both cases are realised [0]. 
In order to proceed, we must allow for the possibility that the average leakage from F to B 
may be different from the average leakage from B to F, and thus introduce different parameters: 
Pf 7^ Vb\ we must also allow for the possibility that, on average, more clans are produced in 
one region than in the other; we call r the average fraction of clans produced in the F region: in 
general 1/2. The symmetric case is recovered for p F = pb and r = 1/2. 

Using the formalism of generating functions, we arrive at formulae for the joint MD in the 
general case. The joint generating function (GF) for each component in Eq. @ is 

G(z F ,z B )= 2J z F F z B B P(n F ,n B ). (6) 

np,riB 

Since clans are Poisson distributed and independent of each other, the joint GF for N F forward 
clans and N B backward clans is 

g{z F , z B ) = exp {N[rz F + (1 - r)z B - 1]} , (7) 
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where N is the overall average number of clans. Calling g c j(z F , z B ) the GF for the joint MD 
of particles from one forward (/ = F) or one backward (/ = B) clan, we make use again of the 
independence among clans and obtain the general formula: 

G(z F , z B ) = G {g CjF {z F , zb),9c,b(zf, z b )) • (8) 

The special case in which the MD is of NB type is of particular interest. In that case in fact 
one has simply (we assume that particles within a clan fall in F or B independently) 

g c ,F(zF, z B ) = g\ og (z F p F + z B q F ) and 9c ,b{.z f , z B ) = g\ og (z F q B + z B p B ), (9) 

where g\ og is the GF for the logarithmic distribution of parameter (3: 

9los (z) =log(l-*/3)/log(l-/3). (10) 

The generating function for the joint probability of n F and n B particles is then given by 

G(z F , z B ) = exp {rN[g Xog (z F p F + z B q F ) - 1]} exp {(1 - r)N[g Xog (z F q B + z B p B ) - 1]} . 

(11) 

This formalism can be used to describe FB multiplicity correlations in very general terms. It is 
interesting to point out a few properties: 

- The marginal distribution is the convolution of two NBD's, thus it can be shown that it is 
an infinitely divisible distribution. 

- When r = 1/2, one obtain for the GF of the marginal distribution a simple correction to 
the NBD generating function; e.g., for the forward marginal GF: 



g F {z) = G(z, 1) 



n i n 2 



1 + ^(1 - Z) + -j^PFQBil ~ Z) 



-k/2 



(12) 



- n F {n B ) is not a linear function of n B , and vice-versa, except for the particular cases 
Pi = 1/2 and pi = 1, (z=F,B); examples are shown in 0]. 

In conclusion, general formulae for the FB multiplicity correlation strength and F and B 
joint distributions generating functions have been given in the framework of the superposition 
mechanism of two weighted MD's for different classes of events. Assuming NB regularity 
behaviour for 2- and 3-jet samples of events in e + e~ annihilation at LEP energy, results obtained 
by OPAL collaborations are correctly reproduced within experimental errors. In pp collisions 
the FB multiplicity correlation strengths in the two substructures (soft and semi-hard events) 
turn out to be quite important: in particular an interesting connection is found between the 
particle populations within clans, particle leakage from clans in one hemisphere to the opposite 
hemisphere and superposition effect between different substructure of the collision. This finding 
favours structures with larger particle populations per clans and the decrease of the average 
number of clans. 

It has been shown that assuming different particle leakage percentages (p B ^ p F ) for bi- 
nomially generated particles from clans in one hemisphere to the opposite one and asymmetric 
(r 7^ 1/2) distribution in the two hemispheres of binomially generated clans, a general formula 
for the generating function of the joint (n F , n B ) -charged particle multiplicity distribution for 
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each class of events (or substructure) in asymmetric reactions can be obtained when the to- 
tal multiplicity distribution is of negative binomial type, and can be easily generalised to any 
discrete infinitely divisible multiplicity distribution. This search is relevant for the study of 
forward-backward multiplicity correlations in non-identical heavy ion collisions and in proton- 
nucleus collisions. Accordingly, the newly introduced particle leakage and asymmetry parame- 
ters can be used for classifying different classes of collisions. 
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